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Abstract 



Let h : [0, oo) [0, cxd) be continuous and nondecreasing, h{t) > if f > 0, 
and m, q be positive real numbers. We investigate the behavior when k ^ oo 
of the fundamental solutions u — Uk of dtu — Au™ + h{t)u'^ = in SI x (0, T) 
satisfying Uk {x, 0) = k5o. The main question is wether the limit is still a solution 
of the above equation with an isolated singularity at (0,0), or a solution of the 
associated ordinary differential equation u' + h{t)u'^ — which blows-up at t = 0. 
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1 Introduction 

Let m and q positive parameters and h : [0, cxd) [0, cxd) a nondecreasing continuous. If one 
consider a reaction-diffusion equation such as 



(m > for simplicity) in a cylindrical domain — x (0,T) {N > 1), the behaviour of u 
is subject to two competing features: the diffusion associated to the partial differential operator, 
here —A, and the absorption which is represented by the term h{t)u''. When q > I and h{t) > 
for t > 0, the absorption term is strong enough in order positive solution to satisfy an universal 
bound 



for every {x,t) S Q^. In addition, the function Uh which appears above is a particular solution of 
(1.1 ). The associated diffusion equation 



dtu - Aw" + /i(t)w« = 



(1.1) 




(1.2) 



dtv - Ai;™ = 







(1.3) 
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admits fundamental solutions v = Vk (fc > 0) which satisfy Vk{x, 0) = k5o if m > {N — 2)^/N. If 

/ / h{t)vl dx dt <oo, Br — {\x\ < R}, (1.4) 

Jo J Br 

for any R G (0, oo], it is shown that (1.1 ) admits fundamental solutions w = Wfe in Q"^ which satisfy 
initial condition Mfc(x,0) = k5Q. The maximum principle holds and therefore the mapping k ^ Uk 
is increasing. If ft, > on (0, oo) then due to universal bound(1.2) there exists Uao = linife^oo U]., 
and solution of (1.1 ) in Q"^. A natural question is whether Ur^, admits a singularity only 

at the origin (0, 0) or at other points too. Actually, in the last case it will imply Uao = U since the 
following alternative occurs: 

(i) either u^a = U. [complete initial blow-up); 

(ii) or Uoo is a solution singular at (0, 0) and such that lim^^o u{x, t) = for all x ^ 0. {single-point 
initial blow-up). 

This phenomenon is observed for the first time by Marcus and Veron. They considered the 
semilinear equation 

dtu - Au + h{t)u'' = (1.5) 

and proved [8, Prop. 5.2] 

Theorem 1.1 If h{t) = e""/* (k> 0), then the complete initial blow-up occurs. 

However they raised the question whether this type of degeneracy of the absorption is sharp 
or not. The method of [8] relies on the construction of subsolutions associated to very singular 
solutions of equations 

dtu - Au + c,t"u'^ = (1.6) 

for suitable a > and > 0, and on the study of asymptotics of these solutions. One the main 
result of present paper states that if the degeneracy of the absorption terms is lightly smaller 
respectivelly to Th. 1.1, then localization occurs. 

Theorem 1.2 If h{t) = ex.p{—u){t)/t), where w is continuous, nondecreasing and satisfies 

/' < 00, (1.7) 

Jo s 

then Uoo has single-point initial blow-up at (0,0). 

The method of the proof is totally different from the one of Marcus and Voron and based upon 

local energy estimates in the spirit of the famous Saint- Venant 's principle (see [5, 12, 13]). Using 

appropriate test functions we prove by induction that the energy of the fundamental solutions 

x 

remains uniformly locally bounded in Q \ {(0, 0)}. 
In the case of equation 

Stu- AM + /i(i)(e" - 1) = (1.8) 
the same type of phenomenon occurs, but at a different scale of degeneracy. We prove the following 

Theorem 1.3 1) If h{t) = e~'^ for some k > 0, then the complete initial blow-up occurs. 

2) Ifh{t) = e-^ ' " for some ui € C(0,oo) positive, nondecreasing and satisfying (1.7 ), then 
has single-point initial blow-up at (0,0). 
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In this paper we also extend the study of equation (1.1 ) to the case m 7^ 1. The situation 
differs completely corresponding to m > 1, the porous media equation with slow diffusion, and 
to {N — '2)+/N < m < 1, the fast diffusion equation. Concerning the porous media equation, we 
prove 

Theorem 1.4 Ifq > m> 1 and h is nondecreasing and satisfies h{t) = 0(f(9-"')/('"-i)) ast ^ 0, 
then Uoo = Uh- 

Wc give two proofs. The first one, valid only in the subscritical case l<m<q<ni + 2/N, is 
based upon the construction of suitable subsolutions, as in the semilinear case. The second one, 
based upon scaling transformations, is valid in all the cases q+1 > 2m > 2 where the Uk exists. 
It reduces to proving that the equation 

-A* - + = in 

admits only one positive solution, the constant 1. The localization counter part is as follows. 

Theorem 1.5 Assume q>m>l,in Equation (1.1 ). Ifh{t) = t(«-™)/(™-i)a;-i(f) withu){t) 
as f — > 0, and 

[ co^{s)—<oo (1.9) 
Jo s 

where 



[Af(m-l) + 2(m+l)](g-l)' 
then Uoo has single-point initial blow-up at 0,0). 

Actually, the method is applicable to a much more general class of equations. 
In the fast diffusion case there is always localization. 



Theorem 1.6 Assume {N — 2)+/N < m <1 and q > 1, in Equation (1.1 ). Then 

Uh{t),cA-^\ ) (1.10) 

where 

l/(l-m) 



\2m{mN + 2- N 



This type of problem has an elliptic counterpart which is initiated in [10] where the following 
question is considered: suppose f2 is a bounded domain in M^, g > 1 and h e C(0, 00) is 
positive. What is the limit, when k ^ co oi the solutions (when they exist) u = Uk of the following 
problem 

f -Au + h(p{x))u'^ = infl 

{ (1-11) 

u = kSo in dil, 

where p{x) = dist (x, 9fi). It is proved in [10] that, if h{t) = e~^/*, then Uoo{-= limfe^oo Uk) is the 
maximal solution of the equation in O, that is the function which satisfies 

f -Au + h(p(x))ui = in n 

\ (1-12) 

[ limp(a.)^o u{x) = 00. 
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On the contrary, if h{t) = t", for a > and l<q<{N + l + a)/{N - 1), it is proved in [11] that 
Woo has an isolated singularity at 0, and vanishes everywhere outside 0. In a forthcoming article we 
shall study this localization of singularity phenomenon for the complete nonlinear elliptic problem, 
replacing the powers by more general functions, and the ordinary Laplacian by the p-Laplacian 
operator. 

Our paper is organized as follows: §1 Introduction. In §2 we study sufficient conditions of complete 
initial blow-up for semilinear heat equation. In §3 we prove sharp sufficient condition of existence 
of single point initial blow-up for heat equation with power nonlinear absorption. In §4 local energy 
method from §3 is adapted to the heat equation with nonpower absorption nonlinearity. §5 deals 
with porous media equation with power nonlinear absorption, §6 — the fast diffusion equation 
with nonlinear absorption. 



2 Complete initial blow-up for semilinear 
heat equation 

We recall the standard result concerning the existence of a fundamental solution u = Uk {k > 0) 
to the following problem 

f dtu -Au + g(x, t,u) = in Q'^ = R^ x (0, T) 

{ (2-1) 
I u{x, 0) = k6o- 



If V is defined in , we denote by g{v) the function {x,t) i— > g{x,t,v{x,t)). By a solution we 
mean a function u € Lj^ciQ^) ^^^^ ^^at g{u) G L]^^{Q^), which verifies 

{-udt(i) - uA(j) + g{u)(j)) dxdt = k(j){0, 0), (2.2) 



for any (f> e Co'^(M^ x [0,T) x R). We denote by E{x,t) = (47rt)--^/2e-l^l'/4t the fundamental 
solution of the heat equation in Q°° , by -Bij(a) an open ball of center a and radius R, and -Bk(O) = 
Bji- The following result is classical 

Theorem 2.1 Let g G C(R^ x [0, T] x M) such that g{x, t,r)>0 on x [0, T] x M+, and assume 
that g = gi + g2 where gi and 52 are respectively nondecreasing and locally Lipschitz continuous 
with respect to the r-variable functions. Let k > be such that 



fl 

Jo Je 



T 

g{x,t,kE{x,t))dxdt < 00. (2.3) 

Br 



for any R> 0. Then there exists a solution u = Uk to problem (2.1). Furthermore, if g2 = 0, then 
Uk is unique. 

Function g(x,t,r) = e~'*/*|r|'~^r, with «; > and q> 1, satisfies (2.3). Thus the problem 

{ (2-4) 

[ u(x, 0) = kSi}. 

admits a unique solution. The next result is proved in [8] , but we recall the proof both for the sake 
of completeness and to present the key-lines of the method in a simple case. 
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(2.7) 



Theorem 2.2 For k > 0, let Uk denote the solution of (2.4 ) in Q°°. Then Uk 1 Us as k ^ oo, 
where 

/ ft N 1/(1-9) 

Us{t) = (^{q - 1) y e-'^^'dsj , Vi > 0. (2.5) 

Proof. Case 1. 1 < g < 1 + 2/N. For any e > 0, u/j = u satisfies 

dtu - Am + e-''/'M« > (2.6) 

on Q^. Therefore \iv = Vk is the solution of 

f dtv - Aw + e-'^/^v'i = in Q°° 
[ v{x,0) = kSo, 

there holds Uk> Vk- Passage to the limit k ^ oo, yields 

lim Uk ■■= Uoo >Voo = lim Vk in Q". (2.8) 

If we write v^{x,t) = e'*/"'(3~^^t~^/(«~^\f (.x/Vi), then / is radial and satisfies 

/"+(^— +2j/' + ^/-/'' = on(0,oo), 
/'(0) = 0, lim,^^r2/9-i)/(r)=0. 

Furthermore the asymptotics of / is given in [2] , 

/(r) = Cr2/(9-i)-^e-'^'/4(l + o(l))) , as r ^ oo, 
for some C = C{N,q) > 0. Therefore 

f{r) > C{r + if/ii-^)-N^-r'/4 > (2.9) 
for some C = C{N, q) > 0. If we take t = e, we derive from (2.8 ) 

Woo(a;,t) > e'^/*(«-^H-i/(«-iV(-jVv^) inK^. (2.10) 



Let < £ < 2^K/{q- 1). Inequalities (2.9 ) and (2.10 ) imply 

«cx>(a;,t)>Ci-V(«-i)e(''/(''-i)-^V4)t-^ Var e S,. (2.11) 

Therefore \imt^QUoo{x,t) = oc' , V.t G Bf . We pick some point xq in Bg. Since for any fc > 0, 
the solution Uks^g of (2.4 ) with initial value kSxg can be approximated by solutions with bounded 
initial data and support in Bcr{xo) (0 < cr < £ — |a;o|), the previous inequality implies 

Woo {x,t) >Uoo{x- XQ,t). 

Reversing the role of and xq yields to 

Uooix,t) =Uoo{x - XQ,t). 

If we iterate this process we derive 

u^{x,t) =u^{x-y,t), Vj/eM^. (2.12) 
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Since Uks^ is radial with respect to y, (2.12 ) implies that Uao{x, t) is independent of x and therefore 
it is solution of 

f z' + e-''/^z'i = Q on(0,oo) 

{ (2.13) 
linit^o z{t) = oo. 

Thus Uoo = Us where Us is defined by (2.5 ). 

Case 2. q>l + 2/N. Let a > such that q < qc,a = 1 + 2(1 + a)/N. We write e"'"/* P/i(t) 
with h{t) = t-"e-^/*. The function h is increasing on (0, K/a] and we extend it by h{0) = 0. Let 
< e < n/a, then the solution u = of (2.4 ) verifies 

dtu -Au + h{e)ru'i > 0, 

/- \ -1/(9-1) 

in X (0,e]. As in Case 1, u is bounded from below onM^ X (0,e] by /i(e) Voo where 

Woo = w is is the very singular solution of 

dtv -Av + t^vi = 0. (2.14) 

Then Vooix,t) = i-(^+")/(«-^)/c«(|a;| /^ft), and /„ = / satisfies 

/"+(^ + 0/' + ^/-/'' = O on (0,00), 

/'(O) = , lim^^oo ^2(l+a)/,-l)y(^) ^ 

The asymptotics of fa is given in [9] 

fair) = c;r-2(i+«)/(«-i)-^e-'''/4(i + o(l)) as r ^ 00, 

thus 

fair) > Cil + r)2(i+«)/(</-i)-We-'''/^ Vr e K+. 

Consequently 

uix, t) > (5e(«/(9-i)-^'/4)t-^ ^ va; e Be. (2.15) 
Taking again < £ < 1\/ iq — 1), we derive 

lirn u(a;, t) = 00 , Mx & B^. 

As in the Case 1, it yields to u^oix, t) = Uooix — y,t) for any y e , and finally Uooix, t) = Usit). 
□ 

Next we consider Cauchy problem for diffusion equation with an exponential type absorption 
term 

f dtU -Au + /i(i)e" =0 in 

\ (2-16) 
[ uix, 0) = k6o 

where h G C(M+) is nonnegative. Theorem 2.1 yields the following existence result: 
Proposition 2.3 Assume h satisfies 

\imt'^/'^lnhit) = -oo. (2.17) 
Then for any k > problem (2.16 ) admits a unique solution u = uu. Furthermore 

Ukix, t) < Vsit) := - In ^ his)ds^ V(a;, t) e Q°°. (2.18) 
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Notice that estimate (2.18 ) is a consequence of the fact that Vs satisfies the associated O.D.E. 

y' + h{t)ey = in(0,oo), 

with infinite initial value. Our main result concerning nonexistence of localized singularities for 
equation (2.16) is 

Theorem 2.4 Let h{t) = e"''"^' for some a > and any t > 0. Then Uk't Vs as k ^ oo. 

Proof. Step 1. Construction of an approximate very singular solution. For n > 1 and c„ > to 
be defined later on, let v = Vn be the very singular solution of 

a*?;- Aw + c„t""w" = 0. (2.19) 

The necessary and sufficient condition for the existence of a is 

n< l + N{an + l)/2. 

This function is obtained in the form 

Vn{x,t) = r(i+"")/("-i)i^(x/v^), 

where F solves 

Ai^ + h.DF + i±^F - c„F" = 0. 
2 n—1 

We fix 

—^ = l + -^ar, = {2 + N){n-l)/2-l, (2.20) 

77- -L ^ 

and set 

fn = cH^-'^F. 

Then /„ solves 

1 + 2 

A/„ + + ^-/« - /" = 0- 

We prove that /„ has an asymptotic expansion essentially independent of n, in the following form 
/40 > + =^ Vn{x,t) > 5c-V(-i)i-2-^/2(|^|2 + f)e-klV4t (2.21) 

It order to see that, we put 

2 ^ 



■^"~ViV + 2j 
then 

1 - Ar + 2~ N + 2 ~ 

^fn + l^^.Dfn + -^fn ' fn = 0- 

By the maximum principle < < 1 so that < /" < /" for n' > n. Thus 

A/„ + l^.Du + ^fn - ^f:' > 0, 

which implies that /„ is a subsolution of the equation for /„< and therefore, 

/TV + 2\ ("'-")/("-!)("'-!) 

n'>n^fn<fn'^fn<[ fn'- (2.22) 
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In the particular case n = n* = {N + A)/{N + 2), the equation faUs into the scoop of Brezis- 
Peletier-Terman study since it can also be written in the form 

and their asymptotic expansion applies (with 2/(n* — 1) — A'' = 2) as |^| — > oo: 

fn'iO = C\i\'e-\^\'/\l + o(l)) =^ /„.(e) > S^l^f + l)e-l«l'/4 (2.23) 
Combining (2.22 ) with n = n* and n' replaced by n, and (2.23 ), we get 

2 \ (n-n*)/{n-l){n* -1) 



/n(0>^*(^^J (ICr + l)e-l^l Ve (2.24) 

Since n i— > {2/{N + 2)^"~" is bounded from below independently of n > n*, we get 

(2.21 ). 

Step 2. Som,e estimates from below for a related problem. In order to have w„ < u in the range of 
value of It, which is 

u{t) < Vs{t) = - In h{s)dsj Mt > 0, (2.25) 

wc need u = u„ to be a subsolution near t = oi the equation that u verifies. Furthermore this 
can be done up to some bounded function. It is sufficient to have 

c„i"-(x" + 1) > /^(^)e^ yt e (0,T„], X e [0,Vs{t)] (2.26) 

where r„ has to be defined. In particular, at the end points of the interval, 

'{i)Cnt'"' >h{t) 

,„),„<.,.(„,.f 1 V A, '.(0 P.27) 



We write (2.26 ) in the form 



and set 



Joa'{s)dsJ J jQh{s)ds 



< (2.28) 



1 + ~ h{t) 
<i>{x) 



Then 



(t)'{x) = 



1 + a 

1 + x" — nx' 



n-l 



(l + a;")2 

The sign of (j)' is the same as the one of tp{x) = 1 + a;*^ — nx"~^, a function which decreasing then 
increasing, is positive near 0, vanishes somewhere between and 1 and again between n — 1 and 
n. The first maximum of (f> is less than e/2. This is not important in (2.28 ) since wc can always 
assume that the minimum of Ckt"'*' /h{t) is larger than e/2. Therefore, it is sufficient to have 

< (2.29) 



1 + V^{t) - h{t) ' 
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in order to have (2.28 ). This is exactly (2.27 )-ii. If we express h{t) in the form 
then (2.27 )-ii is equivalent to 

c„i"" (a;"(<) + 1) > ~Lo'{t). (2.30) 

Since 

a;"(i) + 1 > 2i-"(a;(t) + 1)", 
we associate the following O. D. E. on 



(77 + 1)"' 
the maximal solution of which is 

- i. ( t ] .-(a„ + l)/(n-l) _ i ( ^_ \ .-l-N/2 

2\cn{n-l)J -2[cn{n-l)J 
If we write w in the form 

with a(0) = 00, q' < 0, then (2.27 )-ii becomes 

c„t°"' (e""(*) + l) ^ -a'(i)e"(*\ 

and this inequality is ensured provided 

c„t«"e("-i)"(*) > -a'(t) ^Cn> _a'(t)e(i-")"(')-«"'"* = -to'(t)e(i-"H"W+2"'(^+2)i"0, 

(2.31) 

by replacing a„ by its value. Next we fix 

a{t) = aa{t) = j > (2.32) 

where cr > is a parameter, thus 

-iQ;'(t)e(^""H"(*)+2"'(^+2)li*) = g(l-nWt-(2-i(n-l)(JV+2) + l)lnt ^ gp(t) _ 

In order to have (2.31 ) it is sufficient to have the monotonicity of the function p and 

a{n - 1) n{N + 2)-N 



P'{t) 



t2 2t 



Then there exist 7 > 0, independent of k and a such that p'{t) > on (0,0-7]. Consequently, 
inequality (2.31 ) is ensured on (0, e] C (0, (77] as soon as 

C„ > e''(") = ^(l-n)<j/c-2-\n{N+2)-N)\ne_ ^2.33) 

Step 3. Complete initial blow-up for a related problem. Assume now 

h{t) = at-^e^'"-'"' (2.34) 
for some ct > 0. For n > 2, we fix e < 0-7 and take c„ = e''^*^. On (0, e] we have 

c„t""(e""(*) + 1) > -a'(t)e"(*). 
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Therefore, ii u = Uk is the solution of (2.16 ) with h{t) given by (2.34 ), it satisfies u{t) < Vs{t), 
where Vs is given by (2.25 ), and 

dtu -Au + Cnt"" (li" + 1) > in Q^ 
Therefore u is larger that the solution v = Vk oi 

dtv - Aw + c„t"" (w" + 1) = in Q^ 

with ■Ufc(O) = k5Q. Furthermore Vk > Vk — Cnt°'"'^^ / (ocn + 1), where v = Vk solves 

afU-Au + c„t""i;" = inQS 

with i'fc(O) = USq. If we let /e ^ oo, we derive from (2.21 ) and by replacing c„ = e''^'^) by its precise 
value e(i-")<^A-2-'("(^+2)-JV)ine^ ^hat 

C^t"""'""'" r, ,r;o 9 <r I (n(iV+2)-Nlne |a:|2 

Q;„ + 1 

on (0, e]. In particular 

uUx,e)>5e-^-'"\\x\^ + e)e^+-^^^ ^. (2.35) 

Taking < <j/4 yields to 

lim e"^"^/^( a; + e)e"+ ^t^^t^ = oo. 

Thus 

lim Moo (a;, e) =00, Vx e B r^i^- 
As in the proof of Theorem 2.2, it implies Uoo = ^s- 

Step 4- End of the proof. Since for any a > tr > there exists an interval (0, 9] on which 

any solution of (2.16 ) with h{t) given by (2.34 ) is a subsolution in of the same equation with 
h{t) = e~'^ '^^ . This implies the claim. □ 

3 Single point initial blow-up for semilinear 
heat equation 

We consider the following Cauchy problem 

f dtu- Au + h(t)\u\'^~'^ u = mQ°° 

{ (3-1) 
u{x, 0) = kSo- 

The first result dealing with the localization of the blow-up that we prove is the following. 

Theorem 3.1 Assume h{t) = e"'^'^*^/* where uj € C([0, 00)) is positive, nondecreasing function 
which satisfies (jj{s) > s"" for some ag G [0, 1) and any s > 0, and the following Dini like condition 
holds: 

r^ds< 00. (3.2) 
Jo s 

Then Uk always exists and := \im.k^oo Uk has a point-wise singularity at (0,0). 
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Proof. The proof is based on the study of asymptotic properties as fc ^ cjo of solutions u = Uk of 
the regularized Cauchy problem 

r wt- Au + /i(t)|u|9"iM = inQ^, 

i , (3-3) 

i u{x, 0) = uo,k{x) = Ml^^k-^/^Sk{x) Vx G R^, 

where Sk G C(R^), supp(5fe C < , 5k 5{x) weakly in the sense of measures as A; — > oo 
and {Mk} is some sequence tending to oo as fc ^ oo fast enough so that 

Ml'^k-^''^ ^ oo as fc ^ oo. (3.4) 

Without loss of generality we will suppose that 

ll'^feWllL(RiV) < coA;^ VfceN, Co = const. (3.5) 



Our method of analysis is some variant of the local energy estimates method (also called Saint- 
Venant principle), developed, particulary, in [12,13,15-17] (see also review in [5]). Let introduce 
the families of subdomains 

0(t) = n {|a;| > r} Vr > 0, 
Q'-(T) = f}(T) X (0,r) VrG(0,T), 
Q^(t) = n(r) X (r,r) Vre(0,T). 

Step 1. The local energy framework. We fix arbitrary A: G N and consider solution ti = Uk of 
(3.3), but for convenience we will denote it by u. Firstly we deduce some integral vanishing 
properties of solution u in the family of subdomains Qr ■= x (r, T). Multiplying (3.3 ) by 

u(x,t) exp ( — - — ) and integrating in Q^, we get 

V 1 + T - r / 



T-r 



t-r 



+ 1^ (P,n|2 + /,(t)|nri)exp(-^^ 



r 



dxdt 



+ z 7^ / exp ( — — ^- — ] dxdt 



2-^ [ \u{x,r)f dx + 2-'^ [ \u{x,r)fdx, (3.6) 

Jq(t) JR«\0(r) 



/0(t) JR«\0(r) 

where r > is arbitrary parameter. Using Holder's inequality, it is easy to check that 

2 

/ \u{x,r)\'^dx<CT^^^h{r)-^ i [ \u{x,r)\'^+^h{r) dx] . (3.7) 

Js.^\n{T) \Js.^\n{T) J 

Here and further we will denote by c, Cj different positive constants which do not depend on 
parameters fc, r, r, but the precise value of which may change from one ocurrence to another. Let 
us consider now the energy functions 

h{r)= [ \D:,u\^dxdt, h{r) = ( h{t)\u{x,t)\'i+'^ dxdt, h{r) = ( \u\^ dxdt. (3.8) 
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It is easy to check that 

dhjr) 
dr 



[ h{r)\u{x,r)\''+^ dx> ( h{r)\u{x,r)\''+'^ dx Vr > 0. 



Therefore it follows from (3.6) and (3.7) 



/ |u(x,T)pda; + /i(r) +/2(r) +/3(r) < CT^^&T^/i(r)"^ (-/2(r))'+i +c / \u{x,r)\^ dx 



n(r) 

Vr > 0, Vr : < r < T. (3.9) 



Next we introduce additional energy functions 



f{r,T)^ / \u{x,r)\^dx, Ei(r,T) = / \D^u\^ dxdt, E2{r,T) = / \u\^ dxdt. (3.10) 
Ja(r) JQ'-{t) JQ'-{t) 

Now we deduce some vanishing estimates of these energy functions. Let fi be some nondecreasing 
smooth function defined on (0, oo), ii{t) > for r > (a more precise definition will be fixed later 
on). Then multiplying the equation (3.3) by u{x, t) exp(— /x^(T)i) and integrating in domain Q^{t) 
with T > (remember that suppuo,fe C {|a;| < fc~^}) we deduce easily 



2-V,.,r(T) + J^,.(r) := 2-1 / \u{x,r)\'' e^^{-iJi\T)r) dx+ 

Jq{t) 

/ {\Va:uf + m^(t)|u|^) exp(-y^2(T)t) dxdt 
<Ii{t)-'^ f (|Vxu|^ + M^(T)|u|2)exp(-/x2(r)t)dsdi Vr > fc-^ (3.11) 



Clearly there holds 

dJ^^rij) _ f /. ,2 



./an(T)x(o,r) 



/ (|Vxw|^ + m^(t)|mH exp(-/u2(r)t) rfsrft 

JdQ(T)x(0.r) 

+ / 2/x/x'(r)|upexp(-/x^(T)t)rfxdi 

-2 / /;i/;i'(T)i (|Va,u|^ + M2(r)|up) exp(-/x2(r)f)(ia;(it. 



Since /i'(r) > 0, it follows from (3.11), 



If we suppose 



d f 
--r^tiAT) + 2 / h{t)ii'{t)\u\^ exp{-ii^{T)t) dxdt 
dT ^Q'-(r) 

. 2m'(t) 



. (3.12) 



we derive from (3.12) 

/^,r(r) + J^,r(T) < -2M(r) 



>2-\ (3.13) 
WJ,^,r('r) 



dr 



13 



It is easy to check that this last inequahty is equivalent to 

^ exp (£ ^ ds) UAr) < -I {j,Ar) exp (£ ^ rf.) ) Vr > > fc-. 
By integrating this inequality and using monotonicity of the function Jh^t) we get 



Since 



(/: 








m(t) 




2 





ds I =3— exp 



dT 



it follows from last the relation 



/;i,r(T2) 



exp 



Tl 



ds - 1 



'//n,r(T2)exp 



) < J^,r(Ti) Vr2 > Tl > A; ^ 



Now we have to define ijl(t). Let £ > and 

^(r) =£r-^(T-fc-^) Vt>A;-^ 
One can easily verify that condition (3.13) is equivalent to 

Now from (3.14) follow two inequalities 



(3.14) 
(3.15) 

(3.16) 



A{t2) := 



X exp 



e ((t2 - k-^ f - (n -fc-^)^) ^ £^(r2-fc-i) 



4r 



Vr2 > n > fc-i + 2£-VV/^ (3.17) 



and 

/(r,T2)<A(ri) 



exp 



4r i 



exp 



e'(r2-fc-i)2 



Vt2 > Tl > k-^ + 2e-i/2^i/2_ (3_;l8) 



In particular, for £ = 8 ^ we obtain from (3.17) and (3.18), 



Q-(r) 



|V,,m|' 



{T-k-^f 

64r2 



u\ ) da;df < eexp 



64r 



WxW + —r- dxdt 
2r 



Vt > T^'''(r) := fc-i + 4^2^?, (3.19) 



and 

, e2 / (T-fc-i)2 
/(r,T)<^exp 



+ ^ ) dxdt Vt > ^^\r) := fc"^ + 8\/F, 



(3.20) 
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In order to have an estimate from above of the last factor in the right-hand side of (3.19), (3.20), 
we return to the equation satisfied by u, multiply it by the test function Uk{x,t) ex.p {—t) and 
integrate over the domain = x (0, r). As result of standard computations we obtain, using 
(3.5), 

/ \uk{x,r)fdx+ / {\Va:Ukf + \ukf + hit)\uk\'^+'^)dxdt 

< c ||wo,fc||L2(RiV) < cMk ^ oo as A: ^ 00, yr<T. (3.21) 
Due to (3.20), (3.21) it follows from (3.9) 

/ \u{x,T)\^ dx + h{r) + hir) + hir) 

< cir^^^^fe(r)~iTT(-j^(r-))5TT + caMfcr"^ exp ^- ^ ^ VT>7^''^(r). (3.22) 

Relationships (3.19), (3.20) due to (3.21) yield: 

f{r,T)+E,{r,T)+^^^^^E2{r,T)<C2Mkr-'exp(^-^^^^^^ yT>Tt\r). (3.23) 

Step 2. The first round of computations. Next we construct some sequences {tj}, {rj}, j = 
k,k— 1,...,1. First we explicit the choice of from condition (3.3), let namely 

Mfe = e^'. (3.24) 
Then we choose rfe, such that the following relation is true, 

C2r-iexp(^-^^Mft = M^"«, < eq < e'^ (3.25) 

where C2 is from (3.22), (3.23). As consequence of (3.25) and (3.24) we get 

Tk = 8r^/' [(1 - co)e*= + Inr^i + Inca] . (3.26) 
In inequality (3.22) we fix r = + /c"^, then due to definition (3.25) it follows from (3.22), 

/ \u{x,T)\''dx + h{r)+l2{r)+l3{r) 

, JV(g-l) 2 , 2 

< ci(fc-i +Tfc)^?+^/i(r)"^(-/^(r))^ +Mf Vr:0<r<rk. (3.27) 
Ii{r), l2{r), Isir) are nonincreasing functions which satisfy, due to global a' priori estimate (3.21), 

/i(0) + /2(0) + /3(0) < cMk. (3.28) 

Let us define the number rk by 

rk = sup {r : h{r) + hir) + h{r) > 2Mf } . (3.29) 
Then it follows from (3.27) the following diflFerential inequality 

hir) + hir) + hir) + / \uix,T)\'^ dx < 2ciiTk + k-^)^!^ hir)'^ {-l!^{r))^ 'ir<rk. 

(3.30) 
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Solving it, we get 

h{r) + hir) + hir) < cs{Tk + k-^)"" H{r)-^ Vr < r^, 

where 

rr / 2 \ 2/(9-1) 

H{r) = J h{s) ds and C3 = - — - ) (2ci) 

Next we will use more specific functions 

h{t) = exp 



(3.31) 



(g+i)/(9-i) 



q-lj 



t 

where uj{t) is nondecreasing and satisfies the following technical assumption 
i"" < uj{t) < Wo = const Vi : < t < to, < ao < 1. 
It is easy to show by integration by parts the following relation 



(3.32) 



exp 



1 - Sir) r2 
dt > ^-4 — exp 



(1 - ao)a uj{r) 



au)(r) 



Vr > 0, 



where ^(r) ^ if r — > 0. Therefore 



H(r) > c , s. h(r), c = const > 0. 
uj{r) 



(3.33) 



As a consequence we derive from (3.31), using (3.26), 



N 



2(j{r) 



h{r) + Hr) + h{r) < a [8r| ((1 - eo)e'= + Inr^^ + Inca) ^ + 

2 

Lo(r) 5-1 
X 4 — exp . 

Comparing (3.29) and estimate (3.34) we deduce that rj, satisfies 

Tk < bk, 

where bk is solution of equation 



Vr < rfe. (3.34) 



(3.35) 



bl [{l-e,)e''+\nbl'+\nc,Y+k-'\ u:{bk)^^ b'/-^ exp ((|^^ 



= 2M^° =2exp(eoe'=). 



This equation may be rewritten in the form 
2 



lnc4 + 



q-l 



■In 



bk 



a - 1 bk 



+ N\n 



((1 - £0) exp k + Inb^^ + lnc2) ^ + k-\ 



= ln2 + £oe'= VfceN. 



(3.36) 



16 



Since s ^ Ins ^ as s ^ oo, it follows from equality (3.36) that 

2 Lu{bk) 



(1 + c7(fc))eoe'' > Ak + 



9-1 h 



:= iVln 



86^^^ ((1 - £o)e'= + In + In C2) ^ + k-\ ^ 



+ > (l-7(fc))£oe'= VfceN, (3.37) 

9-1 Ofe 

where < 7(fc) < 1, 7(/;;) — > as A; — > 00. Keeping in mind condition (3.32), we obtain easily 



6^ \Ak\<c{\lnbk\+k) VfceN. 



Due to properties (3.38), it follows from (3.37) 



ce > 



> die" Mk G N, di> 0. 



As a consequence of (3.39), (3.38) we obtain also 



ln6^^<cfc VfceN. 



(3.38) 



(3.39) 



(3.40) 



Now using estimate (3.39) we are able to obtain suitable upper estimate of Tk- Thanks to (3.35), 
(3.39) and (3.40) we deduce from (3.26) 



Tk < cb]!^ exp ( ^ ) < c exp I j- 



uj{bk) 



2 J \di exp k 



Using again estimate (3.39) and the monotonicity of the function w(s), we deduce from the above 
relation 

1 1/2 

, Wo is from (3.32 ). (3.41) 



Tfc < C 



die^ 



Therefore, from inequalities (3.23) and (3.34), definitions (3.25), (3.29) and property (3.35), we 
derive the following estimates 



h{rk) + Hrk) + hirk) < 2M'^° where rk is from (3.35 ), (3.29 ), 
/(r-fe, Tk + fc-i) + Ei{rk,Tk + fc-i) + ^Mrk, Tk + fc"^) < Ml°, 



(3.42) 



(3.43) 



where Tk is from (3.26 ), (3.41 ). Because eo < e ^, it follows from definition (3.24) of sequence 
Mk that 

3Mf <cMfc_i Vfc>A:o(c), (3.44) 

where c > is arbitrary constant. Therefore, adding estimates (3.42) and (3.43), we obtain thanks 
to (3.44) and the fact that Tk ^ rk (which follows from (3.25)), the inequality 



f{rk,Tk + k-^) + J2liirk) + J2^i(''k^^k + k-^) <cMk-i Vfc > fco(c). (3.45) 



i=l 
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Step 3. The second round of computations. Next we introduce the terms rk-\, Tk-i- Firstly we 
come back to inequality (3.14). Fixing here the function 



n{t)=£r-'^{T-k-'^ -Tk) yT>k-^+Tk 
instead of (3.15) and using estimates (3.16)-(3.20), we obtain 



(3.46) 



(T-fc-l-Tfc)^H^ 

64r2 



< eexp 



dxdt 
(r - fc-i - Tfc)2 



64r 



)[ (\VM'+^-$-)dxdt (3-47) 



\/T>4''-^\r) := fc-i + Tfe + 4^/2^/^, 



and 

fir,T) < ^^exp 



(r-fc-i-Tfc)- 
64r 



Vxwl + V- dxdt 



2r 



Vr > ^o^'^"^^ := fc-i + Tfe + S^/F. (3.48) 



The integral term in the right-hand side of (3.47), (3.48) is estimated now by using estimate (3.45) 
obtained in the first round of computation. So, we have 



/ flVxul^ + — ^ dxdt < (2r)-^ 



< c{2r)-^Mk-i 



Vfc > fco(c), Vr > rfc. (3.49) 



Using this estimate we deduce from (3.47) and (3.48) 
(r - Tfc - fc-i)2 



f{r,T)+E,{r,T) + 



64r2 



-E2{r,T) < C2r ^Mfe_iexp 



(T-rfc-fc -i)2 
64r 



Vr >f^''"'^(r). (3.50) 



This estimate is similar to estimate (3.23) from first round. Now we have to deduce the analogue 
of estimate (3.31). For this we return to the starting relation (3.9), where we now estimate last 
term in right-hand side by estimate (3.48), using additionally (3.49). As a result we have 



y^Jiif) < ciT^q+i^'/t(r) 9+1 (-/2(?'))«+' + C2Mfe_ir ^ exp 



(r-Tfc-fc-i)- 
64r 



yr>rk,yT>Tt~^\r), (3.51) 



which is analogous of estimate (3.22) from first round. Next we define the numbers Tk-i and Vk-i 
by inequalities analogous to (3.26) and (3.29), 



C2r^l,Mk-i exp (-^^) = ^^i, < £o < e'^ 
rk-i = sup{r : h{r) + h{r) + h{r) > 2M^^^ J. 



(3.52) 
(3.53) 
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Now combining inequalities (3.30) and (3.44), and using definitions (3.52), (3.53), we obtain the 
following differential inequality 



X]^»W<2ci(rfe_i+rfc + fc-^)^^^/i(r)-iTT(-7^(r))5TT Vr < rfe_i. (3.54) 



Solving this differential inequality, we obtain an estimate similar to (3.31). Using property (3.33) 
we arrive to 

^/,(r)<C4(Tfc_i+Tfc + fc-i)^^^^|^exp(^^|^^^^ Vr<rfc_i. (3.55) 

As in first round we express from (3.52) Tk-i as function Tk-i{rk~i) (the analogue of (3.26)) 

Tfe-i =8rfe{'i[(l-£o)exp(/e-l)+lnrfe_\+lnc2]'/^. (3.56) 

Inserting this expression of Tk-i into (3.55) and then comparing the obtained inequality with 
definition (3.53), we deduce an estimate similar to (3.35), 

r-fe-i < bk-i, (3.57) 

where bk-i is solution of equation 



C4 



((1 - eo) exp(fc - 1) + In + In ca) + Tfe + fc"^] ^ 



X ^^^'^-l^"'" exp f -M^^_iL^ ^ 2M'^°_^ = 2 exp(eo exp(fc - 1)). (3.58) 
Prom (3.50), and due to definition (3.52), it follows 

/(rfe_i, Tfc_i + Tfe + k-^) + ^"^ E2{rk-i,Tk-i + Tfe + /e^^) + Ei (rfc_i, Tfc_i + Tfe + k-^) < M^i^. 

o4rfc— 1 

(3.59) 

Prom (3.55), due to (3.56), (3.57), (3.58), it follows 

/i(rfe_i) + hirk-i) + hirk-i) < 2Ml"_,. (3.60) 

Summing (3.59), (3.60) and using property (3.44), we deduce new global a priori estimate (the 
analogous of (3.45)) which is the main starting information for the next round of computation 

3 2 

f{rk-i,Tk-i+Tk + k-^)+^Ii{rk-i) + ^Ei{rk-i,Tk-i+Tk + k-^) < cMk-2. (3.61) 

i=l i=l 

We are ready now for the next round of computations, introducing the function 

ll(t) = Sr~^{T - k~^ -Tk - Tfe^i) Vt > k~^ + Tfe + Tfe_i 

instead of (3.46) and estimate (3.61) instead of (3.45). We realize j rounds of such computations. 
As result we obtain 

/ I rfe_,-,^rfe_i +fc-M +Y,Ii(.rk-j)+J2^^ ( rfc_„ ^Tfe_( + fc-M < cMfe_,_i, (3.62) 

V 1=0 J i=l i=l \ 1=0 J 
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which was our main aim. 

Step 4- The control ofrk-j,J2i=o ^k-i as j ^ k with arbitrary fc € N. It is clear that rk-j, Tk- 
are defined by the conditions (see (3.52), (3.53)) 



C2r,-4.M,_,exp(-^)=M-,, 



<eo < e 



rk-j = sup{r : h{r) + h{r) + h{r) > 2M'^^'i^}. 
Similarly to (3.56)-(3.58) we deduce that 

1/2 



Tk-j — or 



(l-£o)e^-^+lnr^^+lnc2 

'^k—j — ^k—j 5 



where bk-j satisfies 



(3.63) 
(3.64) 

(3.65) 
(3.66) 



C4 



1/2 

•^k-j ((1 - eo)e'=-^ + In b^^^ + In ca) + ^ Tk-i + fc'^ 



uj{bk-j) 



i=0 

exp f 



JV 



1,9-1 



(q - l)bk- 



2M^^i =2exp(eoe'^-^). (3.67) 



In the first round of computations we have obtained the upper estimate (3.41) for rfe. Let us 
suppose by induction that the following estimate is true 



Tk-i < C 



di exp(fc — i) 



1/2 



Vi < J - 1. 



(3.68) 



We have to prove that estimate (3.68) holds also for i = j. Obviously condition (3.67) is equivalent 
to (see (3.36)) 



lnc4 + ^ In f ^f^) + ^ . + = ln2 + soe^^ 

q-1 V h-j ) 5-1 bk-i 



(3.69) 



where 



J^i^ = iVln 



iv(<i-i)-4 , 1/2 ^ ^ + S ■^fc- 

b-^^ ((1 - e^)e^-^ + ln(6^^ .) + lnc2) + 



i=0 



Because of the induction assumption (3.68) 



j-i i-i 

i=0 i=0 



Wo 



di exp(fc — i) 



1,(<!-1)JV 

"k-j 



< c / ds := cL, 

Jo s 



therefore 



\Ai''>\<c{\lnbk-j\ + {k-j)+\nL). 



(3.70) 
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From (3.69) due to (3.70) we derive easily 



cek-o > '^^^k-j) > ^^gfc-j yj:k-j>ko = ko{L), 

Ok-j 

where fco < oo do not depend on A;. From (3.71) it follows in particular 

ln6fe^<c(fc-j) yj:k-j>ko. 
Thanks to (3.66) and properties (3.71), (3.72), we derive from (3.65), 

1/2 



(3.71) 



(3.72) 



Tk-j < 8bl/_) ((1 - eo)e'^-^ + Inb^lj + Inca) 

'k-j 



< cb^-j exp 



< -^[a;(6fc_,)]i/2 \fj:k-j>ko{L). (3.73) 



Using again estimate (3.71) and monotonicity of ui{s) we deduce from (3.73) 

1/2 



Tk-j < C 



Wo 



die''-^ 



yj:k-j>ko{L). 



(3.74) 



Thus, we have proved by induction estimate (3.68), for arbitrary k—j > ko{L) with r^, r, satisfying 
(3.66), (3.67) and (3.74). 

Step 5. Completion of the proof. We fix now n > ko{L) and take j = k — n in (3.62). This leads to 



fe-r 



fc-r 



/ k„,^rfe_i + fc-M+^/i(r„)+^i;i r„,^rfe_i + fc-M < cM„_i yn>ko{L). (3.75) 



\ 1=0 

Next we have 

k—n 



1=0 



1=0 



Wo 



di exp i 



-,1/2 



< c 



diexp(,i-l) a;(s)l/2 



rfs — > as n — > oo. (3.76) 



Therefore, for arbitrary small 5 > 0, we can find and fix n = n{6) < oo such that from (3.75) 
follows uniform with respect to fc e N a priori estimate, 

sup/ \uk{x,t)\'^dx+ [ [ {\V ^Uk\^ + \uk\'^)dxdt<C = C{S) <oo VfceN. (3.77) 

t>0J|a;|><5 JO J\x\>5 

Since Uk{x, 0) = V > fc-^ V fc € N, it follows from (3.77) that Uoo (a;, 0) = V a; ^ 0, which 
ends the proof. □ 



4 Regional initial blow-up for equation with 
exponential absorption. 

The local energy method we have used in the proof of Theorem 3.1 is based on the sharp interpo- 
lation theorems for functional Sobolev spaces, which are natural tool for the study of solutions of 
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equations with power nonlinear ities. Here we propose the adaptation of mentroned method to the 
equations with nonpower nonUnearities. 
Thus, we consider the Cauchy problem 

J atu-Au + /i(t)(e"-l) = mQ°° 
\^ u{x, 0) = kSo, 

Theorem 4.1 Assume h{t) = e~'^ * ^' where u € C([0,(X))) satisfies the same asumptions as in 
Theorem 3.1. Then solution Uk always exists and Uoo '■= limk^ooUk has a point-wise singularity 
at (0,0). 

Proof. We will consider the family Uk{x, t) of solutions of regularized problems: 
r Ut-Au + h{t){e^ - 1) = in Q"^, 

I (4.2) 
\ u{x,Qi)=uo^k{x)=Ml'^k-^/'^5k{x) Va;eR^, 

where 5k is nonnegative, continuous with compact support in -B/j-i, satisfies estimate (3.5) and 
converges weakly to (5o as fc ^ oo, {Mfe} satisfies condition (3.2). Let us introduce the energy 
functions (we omit index k'muk): 

I^ o{r)= I \V^ufdxdt, 7,(r) = (g!)"M h{t)\u\''+'^ dxdt, h^o{r) = j \u\'^ dxdt. (4.3) 

Multiplying (4.2) by uix, t) exp f — - — ) , integrating in Qr and using equality 

V l + T - r J 



s{e^-l) = Y^ 

we obtain easily 



/i,o(r) +^/,(r) +/3,o(r) < c(g!)2/(''+V^(«-i)/('+i)ft(r)-2/(''+i) (-j;(r))2/(«+^) 

+ c / \u{x,r)fdx Vr > 0, Vr : < r < T, Vg e N. (4.4) 

■Jn(r) 



1=1 



We introduce the additional energy functions 



/(r,r) from (3.10 ), Eio{r,T)= \D^uf dxdt, E2o{r,T)= {ul"^ dxdt. (4.5) 

Jq-{t) ' Jq-(t) 

Instead of (3.21) we derive the following global a priori estimate: 



/ \uk{x,r)fdx+ [ l\V,u\^ + \ukf + h{t)f]^-^^] 



dxdt 



<cMk yr<T. (4.6) 

Using estimate (4.6) instead of (3.21) in a similar way as in the proof of Theorem 3.1, we obtain 
the following inequality, analogous to (3.23), 
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f{r, r) + ^i,o(r, r) + 54^2 ^2,0 (r, t)+ 



<C2A4r-iexp(-^^-^^-i) Vr > f^'^r) = fc-^ + SV^. (4.7) 



Using this estimate we deduce from (4.4) 

E°° 2 iV(t;-l) 2 2 

/((r) +/3,o(r-) < c(g!)<j+iT «+i /i(r) «+i (-/^(r)) ^+1 



+ C2A4r-i exp ( --^^-^^-^ ) Vr > ^^''^(r), Vg G N. (4.6 



Next, we define the numbers Tk, rk- Firstly, set 

rk := sup I r : /i,o(r) + ^i(^) + ^3,o > 2M^° I , < eq < 6"^ 



(4.9) 



Then we fix the sequence {Mk} by (3.24) again and Tk by inequahties (3.25), (3.26). Thanks to 
these definitions we derive the following series of inequalities from relations (4.8) 



;=i 

(4.10) 

Solving these differential inequalities we obtain the estimates 

■00 

/i,o(r)+^/i(r)+/3,o(r) <C3(rfc+/c-i)^(g!)i^i?(r)-i^ Vr < r^, VgsN, (4.11) 

where H{r) is from (3.31). We have now to optimize estimate (4.11) with respect to parameter q. 
By integration by parts, it is easy to check the following inequality 



H(r) > c-Ct exp ( Mr) V r > 0, c > 0. 
uj{r) V r J 

Using Stirling formula q\ ~ y-j and estimate (4.12), we deduce from (4.11) 

00 

h,o{r) + I^^!(r-) + hAr) < C4(t + fc-')^F,(r) Vr < r^, 

where 



(4.12) 



(4.13) 



1=1 



Fq{r) = q^ijj{r) I-'- r -j-i cxp 



2 w{r) 



g — 1 r 

Fixing here the optimal value of the parameter q: 

'w(r) 



exp 



2 /uj{r) 
exp ' 



9-1 



? = 9 



2 exp 



where [a] denotes the enteger part of a, we obtain easily 

/2w(r) 
J^5-<cexp' ^ ' 
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/i,o(r)+f;/i(r)+/3,o(r) <C5(Tfe + A:-i)^exp('^) Vr < r^. (4.14) 



Therefore it follows from (4.13) 

Comparing now definition (4.9) of and estimate (4.14), and using additionally the expression 
(3.26) of Tfc, we obtain 

rk<bk, (4.15) 



where bk is defined by the equation 



C5 



f>^/'((l - £0)6*= + In + lnc2)i/2 + ^-i' 



N 



exp 



[ Mb 
V bk 



= 2Mf = 2exp(eoexpA:), < Sq < e"^ (4.16) 



By an analysis similar to Step 2 in the proof of Theorem 3.1, we obtain estimates (3.37)-(3.40) 
for bk- Then we prove the validity of estimate (3.41) for rfe. As a consequence of estimates (4.7), 
(4.14), thanks to to definitions (3.26), (4.9) of Tk, rk and the previous estimates of rfe, ru, we get 

oo 

/i,o(r) + ^/i(r) + /3,o(r) < 2M'^'> , 
1=1 

^2 

f{rk,Tk + k-^) + ^i,o(r-fe,Tfe + fc-i) + -^E2,o{rk,Tk + k-^)< . 

k 

Summing these inequalities, and using definition of {Mk} and property Tk ^ r^, we obtain an 
analogue of estimate (3.45), namely, 

oo 

/(r-fe, Tfe + fc"^) + /i,o(rfe) + Ii{rk) + h,o{rk) + Eifi{rk, Tk + k~^) + E2,o{rk, Tk + k~^) < cMk-i. 

1=1 

(4.17) 

Using (4.17) as global a priori estimate instead of (4.6) and providing a second round of compu- 
tations similar to (3.46)-(3.57) we derive a second global a priori estimate analogous to (3.61), 

oo 

f{rk-i,Tk-i +Tk + k^^) +/i,o(rfe_i) +^/i(rfe_i) + hflirk-i) 

1=1 

+ Eifl{rk-i,Tk-i + Tk + k~'^) + E2fi{rk-i,Tk-i +rfe + fc"^) < cMfe_2. 

Repeating such rounds j-times we derive a corresponding analogue of relation (3.62). It is easy to 
see that estimate (3.76) for constructed shifts Tk-i remains valid. This fact, similar to what was 
used in the proof of Theorem 3.1, yields to the conclusion. □ 

5 The porous media equation with absorption 

In this section we consider the following problem dealing with fundamental solutions of the porous 
media equation with time dependent absorption, 

( dtu - Adul'^-^u) + h{t)\u\'i-'^u = inQ^ ^ 
I u{x, 0) = kSo- 
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It is standard to assume that /i > is a continuous function and m, q are positive real numbers. 
By a solution we mean a function u € Lj^ciQ^) such that e L]^^{Q'^), hu'^ G ^Lc('5^) and 

y" y" (-wSt^ii - lup-^uA^;. + da;dt = fc(/)(0, 0) (5.2) 

for any cj) e Co^^(M''^' x [0, T)). If /i = and m > (iV - 2)+/Ar this problem admits a solution for 
any k> 0. When m > 1 this solution has the following form 

(r, \ 1/(to— 1) 
2miV ' ^^-^^ 

where 

^= , ^ andCfc=a(m,7V)fc2(™-i)^/^. (5.4) 

A* (m — 1) + 2 

Since is a supersolution for problem (5.1 ), a sufficient condition for existence (and uniqueness) 
of Uk is 

// Bl{x, t)h{t)dxdt <oo. (5.5) 

By the change of variable y — t^^^x this condition is independent of k > and we have 
Proposition 5.1 Assume m > 1, q > 0. If 

f h{t)t^-^i dt < oo, (5.6) 

^0 

then problem (5.1 ) admits a unique positive solution u = Uk- In the particular case where h{t) = 
0{t°') (a>0), the condition is 

Njq - m) - 2 

" > Nim-l)+2- (^-^^ 

We recall that if g > 1 and m> {N — 2)+/iV, any solution of the porous media equation with 
absorption is bounded from above by the maximal solution Uh expressed by 

,t X -1/(9-1) 



{q - 1) y h{s) dsj . (5.8 



Theorem 5.2 Assume q+1 > 2m > 2 and h G C((0, oo)) is nondecreasing, positive and satisfies 
h{t) = 0(t(9~'")/('"~i)) as t ^ 0. Then for any k > Uk exists and limfc_>oo Uk '■= u^o = Uh- 

Proof. We first notice that 

„ „ , . q — m N{q — m) — 2 
q+l>2m>2^q>m>l and > — 7^ -{ -. 

m-1 7V(m-l) + 2 

Step 1. Case q < m + 2/N. In this range of value we know [14] that there exists a nonnegative 
very singular solution v = Voo to 

dtv - Au" + = in Q^, (5.9) 

and Voo = lim/5^ Vk, where the Vk are solutions of the same equation with initial data kS^. Fur- 
thermore, Woo is unique [6], radial with respect to x and has the following form 

V^{x,t) = i-V(9-l)F(|cc| /f(9-™)/2(9-l))^ 
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where F solves 



^ ^ V ^ ' 2{q-l)' q-1 ^ (5.10) 

F'{0) = and lim^^oo ?72/(«-™)F(?7) = 0. 

Actually F has compact support in [0, ^o] for some > 0- Let 7 = (g — m)/{m — 1), then for any 
e > 0, u = Uoo satisfies, for some c > 0, 

dtu - Am™ + ce'^w« > in Q^ 

If we set We{x,t) = a^v^{x,at) with 6 = l/(m - 1)- and a = e-ic-(«-^)/(«-'"\ then 

dtWe - Aw™ + ce'^wl =0 in g^. 

By comparison Moo > w'e in Q^- If we take in particular f = e, it implies 

If < = c~('"^^)/^*^^~^)^o, we derive that Y\va.t^QUao{x,t) = oo, locally uniformly in B^^. This 
implies Uao = Uh- 

Step 2. Case q>m + 2/N. We give an alternative proof valid for all q. We first observe that it is 
sufficient to prove the result when h{t) is replaced by P. If we look for a family of transformations 
u Ti{u) under the form 

Te{u){x,t) = ru{£'^x,it) y{x,t) €Q°°,W>0 

which leaves the equation 

dtu - A|u|"-1m + = (5.12) 

invariant, wo find a ~ {1 + j)/(q — 1) and (3 — {q — m — ^{m — l))/2{q — 1). Due to the value of 
7, we have /3 = 0. Because of uniqueness and the value of the initial mass 

Ti{uk) = ueck V£ > 0, Vfc > =4- r^(woo) = Uoo V£ > 0. (5.13) 

Therefore 

r {x, It) =u^{x,t) y{x,t) GQ°°,yi>o. 

In particular, if we take £ = t~^, 

Uoo{x,t) = t~"Uoo{x, 1) = t~°'(j){x). 
Plugging this decomposition into (5.12 ) yields to 

-ai""- V - t-°"^A(f)'^ + = 0, 

where all the exponents of t coincide since 

m m , m 

am = , — 7 = and a + 1 = . 

m— 1 m— 1 m— 1 

Therefore ^ is a. positive and radial (as the Uk are) solution of 

-acj) - Ac^™ + </)« = in R^. 
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Setting = (f)"^ yields to 



-At/- 3— V)^/" + V'^™ = in M^. (5.14) 

TO — 1 



Clearly ip = ipQ = {m — l)~'"/('3~i) is a solution. By a standard variation of the Keller-Osserman 
estimate, any solution is bounded from above by xpQ. Putting ■(/'(a;) = yl^(a), it is easy to find 
A> Q and a > such that 

-AVi - ViV™ + ^9/™ = in M^, (5.15) 
with < ■)/' < 1. Writting ij) as & solution of an ODE, we derive 

^{r) = Vi(0) + /'^s^-" /'(V;«/™ - V;i/'")a"-^rfs Vr > 0. 
Jo Jo 

If t/js/"* is not constant with value 1, the right-hand side of the above inequality is decreasing with 
respect to r, and the only possible nonnegative limit is 0, by La Salle principle. Thus 

r 2 

for r > To, large enough. If A'' = 2, we set r = Inr, ^'(r) = tp{r) and get 

for r > Inro. The concavity of ^ yields a contradiction. If A'' > 3, we set r = r^~'^/{N — 2) and 
^'(t) = r'^-'^'tp{r). Then * satisfies 

+ C^^(4-JV)/(JV-2)-l/m^l/m < q 

Again the concavity yields a contradiction. In any case we obtain that ^' = 1, or, equivalently 
= Vo and finally, u^o = t-i/f™-^)^^"- □ 

Theorem 5.3 Assume q > m > 1 and h G C((0, oo)) is nondecreasing, positive. If h{t) = 

^(g-m)/(m-l)^-l(^-) mi/l w{t) ^0 ttSt^O, 



I J{s)— < oo, (5.16) 
Jo s 



where 



[N{m-l) + 2{m+l)]{q-iy 
then Uoo '■= limfe_>(x) Uk has a point-wise singularity at (0, 0) 

Proof. The structure of the proof is similar to the one of Theorem 3.1. We study the asymptotic 
behaviour as A; — > oo of solutions u = Uk{x,t) of the regularized Cauchy problem 

Ut-A(|u|'"-iu) + /i(t)|u|«-iu = inQ"^ 

(5.17) 



u 



{x, 0) = uo^kix) = M^+' k -+1 6k{x) X e 
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where 5k is as in Theorem 3.1. Let us rewrite problem (5.17) in the form 

{\v\P-^v)t - Au + h{t)\v\a-^v = 0, in Q'^ 

v = Vk = |ur"^u, p = 1/m, g = q/m (5.18) 
\v{x,0)\P-'^v{x,0) = \voM^~'^vo,k ■■= uo,kix) = Mf^'k-^5k{x). 
Without loss of generality we may suppose 



p+i 



ll'^'fe(a^)ll/^(R«) = / \h{x)\'^ dx<cok'^ VfceN. (5.19) 
Now sequence {Mfe} is such that 

Mj^A;"3^^oo as k^oo. (5.20) 
Step 1. The local energy framework. Consider the following energy functions 

h{T)= [ \V^vfdxdt, h{T)= j h{t)\v\^+^ dxdt, hir) = f \v\p+^ dxdt. (5.21) 

jQr- JQt- JQv 

Analogously to (3.9) we deduce the inequality 

/ |w(a;,T)|P+ida;+/i(r)+72(r)+73(r) < cr^^^i^/i(r)-fT^(-7^(r))lT^+c / \v{x,r)\P+^ dx 

Vr > 0, Vr : < r < T. (5.22) 

This inequality will control the spreading of energy with respect to the r- variable (the time direc- 
tion). As to vanishing property of energy in variable r, we will use the finite speed propagation of 
support property for porous media equation with slow diffusion. In the domain Q^*") (r) we will use 
the energy function i?i(r, r) = \\/xv\'^ dxdt from (3.12). Since supp'y(-,0) = suppi'fc(-,0) = 

suppiiQ.fe = {x : \x\ < k~^}, multiplying equation (5.18) on v{x,t) and integrating in the domain 
Q'^^^t), t > k~^, wc obtain after simple computations (see, for example [1,4]) the following 
differential inequality 

/ \v{x,r)\P+^ dx + Ex{r,T) < cr^+^-'^i-Vu-P) f _ p+i-ci-^iki-p) ^ ^^ ^3) 

Jn(T) \ dT ) 

yr>k-\yr>0 where ^.= ,^/^^-^,)+if+V„ l-^.= ,,, ^ + 

A^(l - p) + 2(p + 1) ' N{l-p) + 2{p+l) 

Solving this inequality and keeping in mind that i^i(r, r) > Vr > 0, Vr > 0, we deduce easily 

v{x,r) = Vx : > /c"^ +cor^"''i£;i(r,fc"^)^^^'w^ := fc-i +cox(r), Vr > 0. (5.24) 

Here the constant cq > depends on the parameters of the problem under consideration, but do 
not on r and k. Analogously to (3.25) we deduce the following global a priori estimate 



Ic 



V,t;|^ + r->ri + /i(t)|t;|'?+i)d:rdt<c||t;o,fc|li+j^(MK). (5.25) 

Thus, due to (5.18)-(5.20), it follows from (5.25) 

Ei{r,0)<cMk Vr>0. (5.26) 
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Next we come back to the inequality (5.22). Due to (5.24) it ensues from (5.22) the inequaUty 

h{r)+l2{r)+l3{r)<c{k-'^+x{r))^^^h{r)-^{-i;2{r))^ Vr > 0. (5.27) 
Remark that due to (5.26) we have 

X(r-) < cir^"^^ . (5.28) 

Step 2. The first round of computations. Now we have to define Tfc, rfe. First we impose the 
relation (i e )(i ) 

Tk > cirl^^^' fe-^^ ci is from (5.28 ). (5.29) 

Then (5.27) yields to 

7(r) := Ii{r) + hir) + hir) < c(fc"^ +rfe)^w^/i(r)"iTT(-7'(r))tTT Vr : < r < rfe. (5.30) 
Solving this differential inequality we get the estimate 

/(r) < /"^^"'^"""Ci Vr : < r < rfc. (5.31) 



{j;his)ds) 



g—p 



Remember that the function /i(s) has the form /i(s) = ^^/(^ p^lo{s) ^, therefore estimate (5.31) 
yields to 

j^,^^ cMr)^^ik^'+rkr vr:0<r<r,. (5.32) 
Thus, as second relation, which defines our pair Tk, rk, we suppose the condition 

C2W(r-fe)^(fc~^ +Tk)^ 



r 



p+i 

k 



<cMk-i, c is from (5.26 ). (5.33) 



Moreover, we will find the pair r^, r^ such that the following property holds 

fc-i + Tfc < 1. (5.34) 
Then the next inequality is a sufficient condition for validity of (5.33): 

p+l _£±1 

C2w(rfc) s-prj. ^ " < cMk-i, c is from (5.26 ), (5.35) 
and we can define r^ by equality 

i-p 

rk := (^) ''^'a;(rfe)I^M;_^. (5.36) 

Now we have to choose the sequence {Mk}. Namely, we set 

Mk:=e'' V/cGN, (5.37) 

and we define Tk, in accordance with assumption (5.29), by 

(i-ei)(i-p) 

Tk=cirlr"'M^ , ci is from (5.28 ). (5.38) 
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Further, due to (5.36) an (5.37), it follows from (5.38), 



N(l-p)+2(p+l) 



ci(^) coinf, (5.39) 



where S = = (g-p)[Af(f-p)+2(p+i)] • From definition (5.36) and because of (5.37) and 

(3.43), there holds 

rfc< (^)'^a;|^exp(-i^(/e-l)) := cgexp ( - 1^ fc), (5.40) 

and rfc ^ as fc ^ oo. Therefore, since u}{s) ^ as s — > 0, it follows from (5.39) that rfe ^ as 
k ^ oo. Consequently we can suppose k so large that condition (5.34) is satisfied. Thus, we have 
pair (r^, r^) for large k 

Step 3. The second round of computations. As a starting global a priori estimate of solution we 
will use now, instead of (5.25), (5.26), the following estimate 

h{rk) = L^^^\V,,vfdxdt<I{rk)<cMk-i, (5.41) 

which follows from (5.32), due to definition (5.33), (5.36) of Vk- Using property (5.24), estimate 
(5.28) and property (5.29), it ensues from (5.41) 

Ei{r,k-^ +Tk) <Ii{r) <h{rk) <cMk-i yr>rk. (5.42) 

Since v{x, Vk) = Va; : |a;| > k~^ + Tfc we deduce similarly to (5.23) 

, -, , (p+l)(l-9i) 

\v{x, Tk + r) 1^+-^ dx + El {rk + T,k~^ + Tk + t) < cr <f+i>-(i-«iKi-»') 

a(r) 

X i - ^Eiivk + T, k-' +rk+ r)) ''+'-''''-"^>''-''> Vr > 0, Vr > 0. (5.43) 



Solving this differential inequality, wc obtain 

v{x,rk+r) = yx :\x\>k-^ +Tk+coXi{r), (5.44) 

where xi(r) := r^~^i£'i(rfe +r, fc"^ +Tk)~ i+p Vr > 0. But (5.42) implies 

(i-ei)(i-p) 

Xi(r) < ciri-^iM^_/+'' . (5.45) 
Now we define Tk-i, Vk-i- In the same way as (5.29) we impose 

Tfe-i > cir^:?W,_/+'' . (5.46) 
Similarly to (5.30)-(5.32) we deduce 

I{r) < '^Mr)^ik-'j^rk+rk.ir Vr : < r < r. + r._, (5.47) 

f 1—p 
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The second relation for defining the pair Tk-i, rk-i is analogous to(5.33) 



< cMk-2, c is from (5.26 ). 



Supposing that 



k ^ + Tfc + Tk-l < 1, 

we can define rk-i by the following analogue of (5.36) 



. in? 
C2 \ p+i 



»-fe + rfe_i := ( - ) ^^^uj{rk + rk-i) — M^_^+' . 



And in accordance with (5.46) let us define r^-i by 



(5.48) 
(5.49) 

(5.50) 
(5.51) 



— w(rfe + rfe_i) 9-p M,. \ '^'M; 



,Y(l-p) + 2(p+l) 



(l-9l)(l-p) 



Due to (5.50) we have 

Tfc-i < ci [(rfe + Tk-iY+^Mlzl] W^^^tW+tj 
< ci 

/ ec2 \ - — i+p , , , c 
= ci(^— j a;(rfe+rfe_i)- 

where S is from (5.39 ). Notice that, duo to (5.47), (5.48), wc have also 

h{rk+rk-i) < I{rk + rk-i) < cMk-2, (5.52) 

and, analogously to (5.42), 

Ei{r, fc-i +Tk+ Tk-i) < h{r) < h{rk + rk-i) < cMk-2 \/r>rk+ Vk-i- (5.53) 

Step 4- Completion of the proof. Estimates (5.52), (5.53) we can use instead of (5.41), (5.42) for 
third round of computations. After j such rounds we deduce that 



Ii(Y.rk-i) < I(^rk-^^ < cMk-j, 

^ i=0 ^ ^ i=0 ^ 

Ei(r,k-^+J2^k-z) < h{r) < h(Y.rk-^ < cMk-j Vr > ^ 



rk- 



where 



(l-8i)(l-p) , I 

^ 1=0 

with the same S as in (5.39), and 

Ea-. = (f)""(E 



rk- 



1=0 



i-p 

g — p _ 1— p 

M~k^t\ 



(5.54) 
(5.55) 

(5.56) 
(5.57) 



31 



Estimates (5.54) will remain true as long as the following analogue of relation (5.49) is valid 

j 

Now we will check this condition. Due to (3.32), it follows from (5.57) 

Yjk-i < (f ) ^c.o^M;_S := CM;_g, = Cexp ( - i^(fc - i - 1)). 
Therefore, from (5.56), it follows 



s 



Thus we have, using in particular the monotonicity of function co{s) 

(l-p)(A:-^-l)^^T5 



g.._.<C.gKCexp(-li^mZl^))] 

- Jk-j-l^ ^ ^ p+1 ))\ 1-p Ja, s 



Ai=Cexp(-^—^k), A2=Cexp\-^-^(k-j -1)]. (5.58) 

p+1 p+1 

Due to condition (5.16) and estimate (5.58) we can find ko G N, which depends on parameters of 
problem under consideration, but does not depend on fc e N, such that 

k—ko 

^rfe_i + fc-i<l VfceN. 

i=0 

At end, our estimates (5.54)-(5.57) are true for all j < k — ko. Therefore the proof of Theorem 5.3 
follows from estimates (5.54)-(5.57), in the same way as Theorem 3.1 from estimates (3.75)-(3.77). 

□ 

6 The fast diffusion equation with absorption 

When (1 — 2/A'')+ <m< 1, it is known that the mere fast diffusion equation 

dtv - Av"" = in Q°° (6.1) 

admits a particular fundamental positive solution with initial data k6o {k > 0) called the Barenblatt 
-Zeld'dovich-Kompaneets solution, expressed by 



-l/(l-m) 

5.(x,i)=t-MC. + A±-^-^) , (6.2) 



{i-m)e \xf Y 

2mN t'^i/N J 



where £ and Ck are given in (5.4 ). The main feature of this expression is that limfe^oo dk = 0, 
therefore 

/ \ l/l-m) 

Um Bk{x, t) = W{x, t):=cA-^\ , (6.3) 



k^oo 
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where 

2m{mN + 2-N) 

This solution has a persisting singularity and is called a razor blade [18]. It has also the property 
that 

liniW(x,t) = Vx^O. 

This phenomenon is at the origin of the work of Chasseigne and Vazquez on extended solutions of 
the fast diffusion equation [3]. Concerning problem (5.1 ), Proposition 5.1 is still valid provided 

TO > (1 + Wc shall denote by u = Uk the solutions of (5.1 ). Furthermore estimate (5.8 

) holds. Combining this with the fact that the are super solutions for the Uk, we derive the 
following 

Theorem 6.1 Assume (1 — 2/N)_^_ < to < 1 and h G C(0, oo) is positive. Assume also that (5.6) 
holds. Then := limfe^oo ^fe has a point-wise singularity at (0, 0) and the following estimate is 

verified 

,2 \ -1/(1-™) / rt \ -1/(9-1)" 



Uoc{x,t) < min J a*-^ [tmN) ' h{s)ds^ 



(6.4) 



Remark. The profile of u^o near (x, t) = (0, 0) is completely unknown. In particular a very chaleng- 
ing question could be to give precise estimates on the quantity min {PF(a;, t), Uh{t)} — Uoo{x, t). 
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